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AflQIHSITION 



by A. J. Viterbi 


I. INTRODUCTION 


Phase-locked or APC loops have found increasing applications in 
recent years as tracking filters, synchronizing devices, and narrow- 
band FM discriminators. Considerable work has been performed to 
determine the noise-squelching properties of the loop when it is 
operating in or near phase lock and is functioning as a linear 
coherent detector (Refs. 1 and 2). However, insufficient consideration 
has been devoted to the non-linear behavior of the loop when it is out 
of lock and in the process of pulling in. Experimental evidence has 
indicated that there is a strong tendency for phase-locked loops to 
achieve lock under most circumstances. However, the analysis which 
has appeared in the literature (Refs. 3, 4, and 5) is limited to the 
acquisition of a constant frequency reference signal with only one 
phase-locked loop filter configuration. 

This work represents an investigation of frequency acquisition 
properties of phase-locked loops for a variety of reference-signal 
behavior and loop configurations. Results are obtained concerning 
the frequency pull-in and tracking behavior (for both constant and 
linearly varying reference frequencies) of the following loop filter 
transfer functions: 


ttiis paper P res * n * s the results of one phase of research 
carried out at the Jet Propulsion Laboratory, California Institute 
of Technology, under Contract No. NASw-6, sponsored by the National 
Aeronautics and Space Administration. 
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1*. K (a wide-band amplifier or resistive pad) 

2. K jjL + a/p j (a single integration) 

3. K[(p+a)/(p + a^(a simple passive RC filter or 
imperfect integrator) 

4. K + a/p + b/p^ j(a double integration) 

The general approach involves solution of the non-linear 
differential equation which describes the system behavior, using 
analog and graphical methods. The terminology adopted in the text 
is the one in general use for closed-loop systems and servomechanisms. 
[However, the application of phase-locked loops to the filtering of 
inoisy signals has generated an alternate terminology in which the loop 
noise-bandwidth is the key parameter (Ref. l). The results will be 
derived in terms of both sets of parameters and are summarized in 
Sec. IX. 


II. BASIC OPERATION 

Figure 1 represents the loop configuration. The frequency of 
the input reference signal, <p s , may be either constant or time varying. 
The frequency of the output of the voltage-controlled oscillator (VCO) 
consists of a constant equal to the center frequency of the oscillator, 
co c , plus a time- varying term proportional to the actuating signal. 

The VCO output is 90 deg out of phase with the input reference so 
that the multiplier output is proportional to the sum of two sinusoids 
of frequencies equal to the sum and difference of the two input 
frequencies. 
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Since the low-pass filter following the multiplier will not pass 
the sum frequency term, the multiplier output, e^* is proportional to 
the sine of the difference phase, 9. The linear transfer function of 
the loop filter shall be termed K 1 F(p) where F( 0 ) =1 and operational 
notation is used throughout. Then 


c l c 2 

©d = *0 sin 9 = ~o sin 


¥ s (t) - u> c t - 


KqKiK^p) 


sin 9 


CD 


Clearly, the proportionality constant, Kq, is one-half the produet of 
the input and return signal-amplitudes. Thus, 


sin 9 = sin 



co c t 


- SM sin 9I 
P T J 


( 2 ). 


where K — Kq K2 shall be defined as the loop gain. Differentiating 
the arguments of both sides of Eq. ( 2 ) yields 


9 + KF(p) sin 9 = <p s - co c 


( 3 ) 


This then is the general equation for phase error of the phase-locked 
loop. At this point the following approximation is generally invoked: 

sin 9 c* 9 for 9 < < 1 rad (4) 


which permits linear analysis to be performed and particularly the 
small-signal response of the system to noisy inputs to be considered. 
To determine the pull-in behavior of the loop, however, the non-linear 
Eq. ( 3 ) oust be solved. 


Page 3 


.-.M 1 li s t. It ■ ir V. 


w e^ - s: : — t 't 


External Publication No. 673 
III. FIRSTS ORDER LOOP WITH CONSTANT FREQUENCY INPUT 

The first and most simple loop which shall be considered is one 
in which the loop filter has no energy storage elements and thus 
produces a constant transfer ratio (see Refs. 3 and 5). The reference 
input shall be a constant frequency sinusoid. Thus Eq. (3) becomes 
the first-order differential equation: 


<p - Q - K sin <p 

where 

Q = <Ps " °&c = constant 

Figure 2 shows <p plotted against If the VCO is initially at its 
center frequency, «(0) = Q and ,(0) = ± nw ( where n is an integer) . 

If the frequency error * is positive, the phase error <p tends to 
incre.se, for negative *, t decreases. Hence, if „ is an even lnteg „ 

and £2 < K, the system will travel along the sinusoidal trajectory of i 
|Fig. 2 until it reaches the T axis at * * sin* 1 (&*) + no. This is a 
stable point; $ can not become negative because <p would then tend to 
decrease and return the system to the <p axis. If n is an odd integer, 

the system will go through a larger part of the sinusoidal trajectory 
until it reaches a stable point at 

(n + l)o + sin'i^Sj 

If 0 >K, however, the trajectory never crosses the <p axis and a 


, i 
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. stable or phase-lock point is never reached. Thus the pull-in 

frequency range of the loop is = K rad/sec. For the first-order 

loop K is also equal to the closed-loop 3 db bandwidth. 

When phase-lock does occur, the steady-state phase error is 
shown to be sin The linear analysis which follows from the 

approximation (4) arrives at a steady-state error Q/K, which is 
clearly an accurate estimate when Q/K < < 1. 

Pull-in time may also be determined from Eq. (5). Since 
<p = dtp/dt, this can be rewritten as 

dt = 1 

d? Q - K sin 9 

.1 

j whence the pull-in time is given by 


t 


/ 


^ final 


Q - K sin 9 


(6) 


^initial 

’'final cannot taken as sin" 1 (Q/K) + nr since for this value the 
demoninator becomes zero and the integral becomes infinite. This is 


correct since, in fact, the time required to reach the steady state 
is infinite. To get a meaningful result can be taken as 

slightly less than the steady-state phase and the time required to 
reach it can be evaluated. The integral is expressed in terms of 
elementary functions in the standard tables. 

The conclusions to be drawn regarding the first-order loop are 
that it will acquire lock within one cycle provided the initial 
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frequency error is less than the 3 db loop bandwidth and that the 
lock-in time is governed by the initial phase error. 

IV. SECOND- ORDER LOOP WITH CONSTANT FREQUENCY INPUT 

To overcome the limited pull-in range of the first order loop, 
and to provide for tracking linearly varying frequencies with narrow 
bandwidths, an integrator is generally used in the loop filter. Then 

F(p) * 1 + | (7) 

and Eq. (3) becomes 


Kp + aK , . x ^ , % 

pp + ( sin ? ) = q ( 8 ) 

o 

Letting aK * and K = 2( and differentiating, the second-order 

differential equation whleh results is 

t 

p^g + 2£ «OnP (sin g) + co^ 2 sin g * 0 (9) 

| • 

| The characteristic polynominal when the equation is linearized by 

| approximation (4) becomes 

I 

i 

p 2 ♦ 2£ coup + co n 2 = 0 (10) 

In servo terminology o> n is the undamped natural frequency and £ 
the damping factor for the linearized servo loop. 
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Rewriting Eq. (9) without operational notation: 


2 

+ 2C 03 n cos <p ^ + Q3 n 2 sin (p = 0 


Making the substitution 


( 11 ) 


t 


2C co n 


( 12 ) 



oh: 


2 

+ con 


sin 


JX 0 


Normalizing and letting 


(13) 


= « = 1 — 
dx * 2C oon dt 


(14) 


<p + <p cos <p + sin <p = 0 


(15) 


Thus one of the parameters has been eliminated by normalizing the 
time parameter. It would now be of interest to construct a plot 
of ^ as a function of <p as was done in the previous Section. To this 
aim let y - x (the abscissa) and <p = y (the ordinate). This requires 
that y = x and if * y. Then Eq. (15) becomes 


or 


y + y cos 


x + S±S5* 


0 
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y - i ’ £ “ - cos x - ^ (16) 

which relates the slope of the system trajectories to the instantaneous 
frequency and phase errors and permits their plotting. Certain 
observations can be made based on Eq. ( 16) : 

1. For large y (which corresponds to large frequency 
error) the second term on the right becomes small and 
the trajectories become nearly sinusoidal. 

2. For x * 0 (<p * 0), the slope is always - 1. Hence the 
y axis is the - 1 isocline (line of constant slope). 

The overall pattern can be determined by finding other 
isoclines which in general will not be straight lines. 

3. The equation is periodic in x with period 2 it; that is, 
the slope of the trajectories is the same at (x 0 , y 0 ) 
as it is at (xq + 2n7r, Yq) . Hence, to describe the 
behavior it is only necessary to plot the trajectories 
for - rr < x < TT. 

4. At y s 0, x s ± ror (n an integer or zero) the slope 
dy/dx is indeterminate since the second term on the 
right becomes 0/0. These points are called singu- 
larities and are either stable points of the system 
or centers of instability as will be discussed more 
fully below. 

The graphical procedure just .described for finding solutions of 
the non-linear differential Eq. (15) is the so-called "phase-plane" 
method (Refs. 6 and 7). The procedure of determining the isoclines 
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and sketching the various trajectories is a long and tedious process, 
however. On the other hand, Eq. (15) is easily programmed for an 
analog computer as is shown in the simplified diagram of Fig. 3. The 
voltage representing cp is used to actuate the x axis of a plotter and 
that representing $ goes to the y axis. 

This results in the plots of Figs. 4 through 7 for various values 
of C* The 9 = x axis extends from - ir to + tt. The following obser- 
vations may be made from the figures: 

1. For large values of y (positive and negative) the 

trajectories are practically sinusoids. The trajectories 
are traversed from left to right in the upper half 
plane and vice versa in the lower. For large positive 
y a small decay may be noted as the trajectory traverses 
the strip from cp = - ir and 9 = ir. The rate of decay 
becomes greater as y decreases. Since the behavior 
is periodic in each strip of width, 2 ir, the decay will 
continue and increase in each successive strip until at 
x * 9 = n*r (where n is some odd integer) the value of 
y will be below the line A - A. At this point the 
system will stop skipping cycles and the phase and 
frequency error will decay in toward cp * 0, 9 = (n + l)ir 
and phase-lock will have been achieved. The same 
behavior ensues when the initial frequency error is 
negative except that the movement of the trajectories 
is from right to left. The system may be said to be in 
frequency-lock when the errors lie within the lines 
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A - A in a particular strip since when this condition 
exists the loop stops skipping cycles. 

The singular points at x = ±mr (n = 0, 1, 2, ...) are 
indeed points of stability or centers of instability 
as mentioned above. When n is even there are stable 
points where the loop achieves phase-lock. These are 
called nodes or foci depending upon whether the system 
is overdamped (Fig. 4) or underdamped (Figs. 6 and 7). 
At odd multiples of 7 r there is a center of instability 
known as a saddle point. Near this point, no matter 
what the direction of the trajectories maybe, the 
magnitude of the frequency error will decrease until 
the saddle point is almost reached and then quickly 
increase again. Further properties of the singular 
points may be found in Refs. 6 and 7 on phase-plane 
techniques, but these are not necessary for the 
discussion. 

It can be seen by comparing the value of y for which 
the line A - A crosses the point x = - ir on the 
various figures that the region of frequency- lock is 
extended as £ decreases. This, however, is a con- 
sequence of the normalization of Eq. (12). Letting 
the true radian frequency dcp/dt = oj rad/sec, and using 
Eq. (14) there results 

<0 * 2C m n <p = 2C co n y (17) 

Figure 8 is a plot of the limiting value y and cft/ton 
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for frequency-lock as a function of C obtained using 
the analog computer. Clearly, the more damped the 
system the wider is the region of frequency- lock for 
a given loop natural frequency. 

The qualitative intuition obtained by examining the phase-plane 
plots may be put on a somewhat more quantitative basis. First, it 
will be shown that the pull-in range is theoretically infinite. Then 
approximate expressions will be derived for the frequency error decay 
per cycle and the pull-in time. 

If both sides of Eq. (16) are multiplied by y and integrated 
between the limits - ir and + w the result is 




cos x dx - 



vr 


sin x dx 


/- TT 


(18) 


. The second term on the right is zero while the first may be 
integrated by parts to yield: 


1 

2 





sin x dy 


(19) 


If the expression for dy from Eq. (16) is substituted into Eq. (19) 
there results: 


1 

2 


y 2 (n) - y^( -7T ) 




sin x cos x dx 


i 

4C 2 J-TT 


-S- ir P x dx (20) 
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The first term is clearly zero and the second can be rewritten as 


1 

2 


y 2 (ir) - y 2 (-ir) 


7T 


1 - cos 2x 


dx 


( 21 ) 




-7T 


Clearly, for positive y the integrand is positive for all values of x 
between the limits ±rr, making the right side of Eq. (21) negative. 
Therefore, during each cycle of width 2 w, the value of y must decrease 
for all initial values of y, which implies that the pull-in range is 
infinite when a perfect integrator is used in the loop filter. 

To determine approximately the decrease in y per cycle of x it 
may be assumed that for y above the line A - A the rate of decrease 
per cycle is constant. This approximation is, of course, best for 
large y. Then, 


y(x) = y(-ir) - sin x - p x 


(22) 


for - 7 T < x < v where p = 6y/2ir and by is the decrease per cycle, 
Substitution into Eq. (16) yields 


dx 


= - COS X - 


sin x 


4C 


y(-tr) 1 - 


sin x £ 


y(-Tr) (y( 


X 

-irl 


If y(-tr) > > 1 


- 6y = y(ir) - y(-rr)^ — x 

4£ 2 


± f W (sin x + ^ ? « sl " * ) dx 

y(-ir) ' y(-^) y(-ir) ' 
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Substituting Eq. (22) for y, assuming y(«ir)> > 1 and p < <1, and 
preserving only the first and second order terms of the expansion 
yields: *. 


6 1 - 


2C tD n y( -tt) 


f w 

, t sin x 

+ sin 2 X 

L 1 

y(-ir) 

u 

y^(-7r)_ 


dx 




jr 2 

2 C co n y(-Tr) 



(26) 


If 6y is small (which is certainly the case for y >3) the ratio 
6t/6y is a good estimate of the derivative of decay time with respect 
to the magnitude of frequency error. Dividing Eq. (26) by Eq. (23) 
yields : 


dt ~2£ 
dy co n 



(27) 


Integrating between the limits y 0 , the initial error, and y A> the 
value at the line A - A, 

Vo 2 - YA 2 - (t 2 ) + -^2 • -^2 12 

' y A' 2y 0 2 2y A 2 

is an approximation of the time required for frequency-lock when the 
initial frequency error is Q Q * 2£ o) n y 0 rad/sec. Of course, the 
derivative approximation becomes increasingly crude as y approaches 
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Evaluation of the integral yields 


by 


tu ♦ m 

4C 2 y 2 (-») 


Substituting for p from Eq. (22): 


by 






(23) 


This expression checks rather well with analog results down to 
y(-ir) = 3. 

By similar methods an expression may be obtained for the time 
required to achieve frequency- lock (defined as the point at which 
the trajectory dips below the line A - A) . Since 



_ dop _ dx 
~ dr dr 



using Eq. (12) the time per cycle is: 


bt = 


br 

can 


1 

2C con 



(24) 


(25) 
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y.A* However, if y 0 was large to begin with, terms of Eq. (28) are 

^mall with respect to the first, and asymptotically the pull-in time 
will be: 


t Oi 


2C y 0 ‘ 


CO 


n 


Q 0 2 
2£ co n v 


seconds 


(29) 


It should be noted that if the integrator in the loop filter is 
not perfect the pull-in range of the loop is no longer infinite, but 


is proportional to the finite time constant of the filter. This will 
be shown in a later section. 


V. SECOND- ORDER LOOP WITH LINEARLY VARYING FREQUENCY INPUT 

When a phase-locked loop is used as a tracking filter or FM 
discriminator it has to follow an input reference with variable 
frequency. If this variation is due to a doppler shift it will be 
almost linear over a large part of the tracking period. If it is due 
to frequency modulation the maximum slope of the modulation will be 
limitsd by its bandwidth. It follows that a loop which can follow 
a frequency which varies linearly at this maximum rate can certainly 
follow the modulation. 

The input frequency to be considered is then cd s + Dt. For a 
second-order loop the general Eq. (3) becomes: 


P9 + aK (sin <p) = Q + Dt 

P 


(30) 


Page 15 



<5$ W* 


External Publication No. 673 


where Q is the initial frequency error which may be zero if the loop 
was originally in lock when the linear variation Dt of the input began. 
Since the right-hand side of Eq. (30) is the difference between the 
reference and VCO frequencies, the time-varying frequency may just as 
well be that of the VCO, corresponding to a constant drift of the 
oscillator. Following precisely the development and substitution of 
Eq. (8) through (15) of the previous section there results: 


? + 9 cos q> + -^2 sin <p * 2 ^" 2 


(31) 


Again let qp = x and <p = yj the slope of the trajectories becomes 

D ' 


- sin x 


dy 

d^ 


,a> 


* - cos X + 


n 


4C 2 Y 


(32) 


The singular points at which the slope becomes indeterminate are now: 


y = 0, x = sin -1 ± 2 mr 

W / 


and 


y = 0, x = n - sin " 1 ( —^5) ± 2 

v co n ^/ 


nir 


(n — 0, 1, 2, •••) 


(33) 


It is interesting to note that if D, the slope of the excitation, 
equals the square of the natural frequency of the loop, the two sets 
of singularities come together; and if 

2 


D > 00 


n 


(34) 


no singular points exist. 
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It will be shown that inequality (34) is the condition for 
instability. That is, if the doppler or modulation slope ever 
satisfies this condition not only can the loop never achieve lock, 
but even if it was initially in lock it will immediately fall out of 
lock when the excitation is applied. 

The phase-plane plot can be obtained by analog methods using the 
computer configuration of Fig. 3 with the addition of a constant bias 
D/a> n ^ on the first integrator. Figures 9 through 13 were obtained by 
thus mechanizing Eq. (31) for a damping factor ^ = 0.707 and various 
levels of the ratio D/oo n 2 . 

The following comments may be made concerning these graphs: 

1. The singular points are of the same nature with the 
linear excitation as without it. The only difference 
is that the stable point is displaced to the right by 
sin -1 ^D/co n 2 j while the saddle point has moved to the 
left by this amount, as was predicted. Since the 
stable point represents phase-lock, it is seen that 
the steady phase-error of a second-order loop with 
linear excitation is 

E ss = sin’ 1 

The linear approximation discussed in Sec. II predicts 
a steady-state error of D/co n 2 which is valid only for 
small values of the ratio. 

2. The trajectories are decidedly asymmetric about the 

x axis. For low values of D/a> n ^ (Fig. 9), in the upper 
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half plane all trajectories above the line B-B are 
divergent , slowly at first and more noticeably as the 
error frequency increases. On the other hand, all 
trajectories below B-B converge toward a stable 
point either in the same strip or some other strip to 
the right or left depending on whether the trajectory 
is in the upper or lower half plane. The line B-B 
exhibits a periodic behavior (neither locking in nor 
diverging) which is, however, unstable since the 
slightest disturbance will tend to push the system 
either into lock or into complete instability. 

Figure 9 verifies the intuitive observation that 
in order to acquire a positive doppler-shifted 
frequency it is best to make the VCO frequency lead 
the input reference (this corresponds to a negative 
frequency error) and let the increasing doppler shift 
decrease the error until lock is achieved. 

3. As the ratio D/a> n 2 is increased the line B-B which 

establishes the limit between stability and instability 
in the upper half plane descends until at D/oo n 2 = 1/2 
(Fig. 10) it practically coincides with the saddle 
; point asymptotes A-A. Thus, when the ratio reaches 
this value all trajectories in the lower half plane 
are drawn into lock, while in the upper half plane lock 
can never be achieved, unless achieved within the 
original strip of width 2 ir. 
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For yet larger values of D/con^ (Fig. 11) another 
phenomenon becomes apparent. While most trajectories 
which start in the lower half plane still achieve lock 
in some strip, there are some that do not. In 
particular, any trajectory that passes within the 
corridor whose lines of demarcation are the extensions 

i 

iof the saddle point asymptotes A-A will be whipped 
| past the x axis into the upper half plane. The width 
of this corridor in the lower half plane is increased 
as D/co n ^ increases (Fig. 12) and many more trajectories 
become unstable. As the ratio becomes almost one 
(Fig. 13) the corridor becomes so wide that only 
trajectories which began in the immediate vicinity of 
the stable singularity ever achieve lock. 

Finally, when D/co^ £1 the singularities 
disappear and no stable point exists. Thus a loop 
which was originally in lock before the excitation was 
applied will immediately lose lock. The conclusion 
which may be drawn from these results is that a 
second-order loop can acquire a noise-free doppler- 
shifted or modulated signal with certainty only if the 
ratio of frequency slope to the square of the natural 
frequency is less than 1/2, and it can track the signal 
once lock is achieved only until the ratio becomes 
unity. 
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It .is_of inte_i^3t£j^_derlY^ an eppCTxiteate expxessioii lQr the 
periodic iimiting behavior between stable, and divergent trajectories 
for small D/co n (line B-B in Fig. 9). If y is to be periodic 
;y(ir) = y(-ir). 

i 

! Multiplication of Eq. (32) by y and integration between the 
limits ±n yields: 


0 = 


1 

2 


y^(ir) - y 2 (-ir) 



cos x dx + 


2C 



Integrating by parts and making use of Eq. (32): 


0 - I sin x f- cos x + 2 ” 2 ' “?sf) dX + ^2^ "2 


-TT 


4C y co n KyJ 2C °>n 


For small D and consequently large y ay of the periodic behavior this 
leads to the approximation 


-tr 


irD 


< Yav ^ 


co, 


n 


(») 


or 


CD 


n 


Vf rv ? “ 

Yav “ 2D 


This can be taken as an approximate expression for pull-in range, 
terms of the condition becomes 



In 


(36) 


1 V 


Page. 20. 


! , A r , , 

\jrf r ,r opui *\o. • i.i;r.cr.:7'r\ 


*?4 

! ’.1 "• _■ } i - ' ' e^f rt 'A r - * ] 

JBIM. fif D/g> A 2 ‘ a/4i £ *3*7 07. this predicts 

i . 1 £LXo 1 1 X O . ^ t * ^ ^ 1 * "* ; 

i Yqy -^ 2.0. it will be noted that this is a fair approx imation. 

2 

I The approximation improves as D/o^ decreases but it is very poor as 
I 2 

iD/o> n approaches 1/2. 


VI. SECOND- ORDER LOOP WITH IMPERFECT INTEGRATOR 

If a passive RC loop filter 2 is used or if the integrator in a 
second-order loop is not perfect, the loop filter has a transfer 
function of the form: 

F( p ) = 2 JL a (37) 

p + a ' 7 

Substitution of this function into Eq. (3) yields 

3t + K (f”H) sin 9 " is - “c 

Assuming the input frequency to be constant and letting Q * <p s 
this may be written as 

~2 + Jt + cos <P ♦ sin <p * aQ (39) 

at* 

To obtain the form of Eq. (9) which represents a second-order loop 
With perfect integrator, again let K = 2£ o> n , aK = co n 2 . 

2 — 

This was also investigated by Gruen (Ref. 3). Preston and 
Tellier (Ref. 4) considered the simpler filter F(p) = a/(p + o) , which 
is a sub-class of this case. 


(38) 
‘ “c 
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Then 


* dt (a + ^ “n cos *) 


+ co n z sin 9 = afl 


(40) 


If the time is again normalized, t = t/( 2C ton), and the result is 


9 


,2t 


+ COS 


9} $ + “^o sin 9 

/ K 2 


J&. 




(41) 


where 


9 


= 4$. 

dT 


)2C “n; 


SI 


With the convention 9 * y, 9 * x 


f _aeV_ 


= . 

dx l a: <°n 


sin x 


+ COS X + 




n 


4C 2 y 


(42) 


This indicates that the singularities are at y * 0, x * sin* 1 ^eft/a^ 2 j 
and x * rr - sin- 1 (ofl/cDn 2 ) which implies that for 


0 ) * 
Q > -Jl 


(43) 


there is no stable point and hence the system can not achieve lock. 
Condition (43) is the upper bound for phase lock; that is, if the 
input frequency is more than co n 2 /a radians away from the VCO center, 
frequency lock-on is impossible on the basis of the singular point 
behavior. 

However, for lesser Q pull-in may not occur even though a 
stable point exists in the phase-plane. Examination of Figs. 14-18, 
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the analog computer solutions of Eg. 41* boars this out. In all eases 
C * 0.707 and a/2£ = CXJ.. The behavior varies with the value of Q. 

Condition (43) requires aQ/to^ < 1 for lock to occur. In Fig. 14 
aO/<%^ =(X4 and pull-in is quite apparent. It should be noted that 
given o/2C ©„ * 0.1 and C s 0*707 =* 1/ then * 2 */T. Thus, if 

the ¥00 was initially at its center frequency, then the initial value j 
of y * <p 3 Q/2£ oo^j ° 2. On the other hand, if initially the VCO was 

not at its center frequency because it had been tracking a signal of 
another frequency, then the initial value of y is not restricted. It 
seems apparent from Fig. 14 that pull-in occurs for all initial values 
of y when Q/a> n * 2 V~2, 

Figure 16 represents the solution for aQ/con 2 = 0*7 or Q/co n * 

7 V5/2. A limit cycle exists about y av = 2.5 toward which all 
higher trajectories converge. If the VCO is at cd c initially then 
initially y - Q/2£ £% * 3.5 and the system converges to the limit 
cycle and exhibits a periodically varying frequency error as long as 
the input frequency remains constant. If the initial frequency of the 
VCO is not co c , one of three situations may exist (see Fig. 19). 

1. The initial VCO frequency (coq) lies further from the 
signal frequency (©$) than does the VCO center 
frequency (««)• Then the initial value of y > Q/2£ a> n 
and the trajectory will converge to the limit cycle. 

2! ojq lies nearer to co s than does co c ; then initially 
0 < y< £1/2 C and the system will lock in if y is 

sufficiently small, or otherwise converge to the limit 
cycle. Figure 16, in fact, shows this behavior in fora 
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of an unstable limit cycle below the aforementioned 
stable one. If y initially lies below this unstable 
limit cycle, lock-in occurs. If it lies between the 
limit cycles it will converge to the higher one. 

(3) coq lies on the other side of co s from co c ; then initially 
y is negative and the system seems to lock in always. 

Figure 15 (aQ/co n 2 =0.6 and £/co n = 3 V~2) represents a behavior 

intermediate to those already considered. Pull-in occurs throughout 

but the rate is slower near the region in which a limit cycle occurs 

for the case of Fig. 16. Thus, it appears that if the VCO is initially 

at co , for C = 0.707 and a/2£ co n = OJ. the limit for pull-in lies in the 
c 

region 

3 ' / ‘ 2< \ Si ( 44 ) 


Figures 17 and 18 represent the solution for Q/co n = 4 J~2 and 
9/2 2. If the VCO is initially at co c , the initial values of y are at 

4 and 4.5 respectively and pull-in does not occur, of course. The 
stable limit cycle is again in evidence and furthermore, if y is 
initially negative (case 3 above), lock-in is not assured since those 
trajectories which pass through a strip determined by the asymptotes 
of the saddle point will converge to the limit cycle. 

An approximate analytical derivation of the pull-in range follows. 
Returning to Eq. (40), letting the frequency error dcp/dt = co and 
dividing by co yields: 


dco _ 

cfep 


(a + 2£ co n cos <p) + 


aQ - co n 2 sin cp 


co 


(45) 
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Multiplication by a> followed by integration from <p = - 7 r to + 7 r yields 
measure of the decay or rise in frequency error during one cycle. 


2 




f° 7 r 

2naQ - f co(cc + 2 £ ao n cos 9 ) dcp 

J’* 


(46) 


Integrating by parts using Eq. (45) 


1 

2 


CD 


^( 7 t) - co^(-tt) 


= 2iraQ - na 


Co( 7T ) + co(-ir) 


(47) 


f (a«f> + 2£ co n sin cp) ^-a - 2£ co n cos <p + 

7 r ^ 


dcp 


Making the assumption that the frequency error co is large compared 
to co n , the terms in go in the integrand will not change appreciably 
during one cycle, so that throughout the cycle co - go(-tt) for -tt< cp < n 
Then: 


1 

2 


<o 2 M 


CD- 


■ 2 ( -*) 


2tt s aQ - 



2 


co 


n 


(a + Cto n ) 
co(-tt) 





(48) 


The trajectories will always decay when the right hand side of this 
equation is negative, which occurs for the following range of Q: 3 


Gruen’s result for pull-in 
technique for co^/a > > 1. It 

(49) by a factor of V^2. 


range was obtained by an empirical 
is more pessimistic than condition 
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2 < 2 “n V 2 + 1 ( 

For the oases investigated «/2S = 0.1 so that the approximate 

condition for pull-in when the VCO is initially at its center 
frequency will be 

^ < 2^6 = (3.46) ^2 

This lies in the region determined on the analog computer (Eq. 44). 

VII. THIRD- ORDER LOOP TRACKING BEHAVIOR 

Section V has demonstrated the limitations of the second-order 
loop in tracking frequency- variable signals. It is of interest to 
investigate the possibility of extending some of these limitations 
by inserting a second integrator in the loop filter. It develops 
that pull-in behavior is less stable for a third-order loop than tor 
one of second order but it will be shown that the tracking range tor 
a loop initially in lock can be extended by means of the second 
integrator. The loop filter transfer function is 

F( P>= 1+ f + ^| <5C 

and Eq. (3) becomes 

o 

n® + Kp ± aKp ± h K / . v 

^ ^2 (sin <p) - Q + Dt (51 
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Letting K = 2£ co n ; aJC = co n 2 and differentiating 

p^qp + f 2C co n p + co n 2 + — b ^ ( sin cp) = D (52) 

Making the substitution 

t _ — 2 — and dropping the operational notation 
2 C co n 

K 2 » n 2 + 4 C 2 “n 2 cos 9 (j*) + “n 2 sln 9 (53) 

+ b 

Normalizing and letting 4jL - 9 

0+9 cos 9 + “3 sin 9 + — 2 / sin 9 ^ = “3 5 

w Y 4£ 2 co n ^ J K un 

Conventional phase-plane techniques are not applicable to third-order 
differential equations because there are three initial conditions 
corresponding to each of the three dynamic variables to be considered: 
phase, frequency, and frequency rate (or displacement, velocity, and 
acceleration in mechanical systems). In principle, the trajectories 
of a third-order equation could be described in three dimensions. 
However, an attempt to project these trajectories for a multitude of 
initial conditions on a plane will generate such a confused diagram 
that little of a general nature can be deduced. 


J sin 9 dx = D 
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On the other hand, if one limits oneself to a single set of 
initial conditions, a meaningful projection of the trajectory onto the 
9 - <p plane will result. The set of initial conditions of greatest 
interest is <p = dep/dt = 0$ d 2 <p/dt 2 = Dj that is, the loop is initially 
in lock so that the frequency and phase errors are zero when the 
reference frequency begins to change linearly. 

The analog computer program of. Fig. 3 can easily be modified to 
include another integrator. Figure 20 represents a series of 
trajectories projected on the <p - <p plane. In all cases C = 0.707 
and D/co n 2 * l/2. In the hypothetical three dimensional "phase space" 
jthe trajectories initiate at $ * D/4C 2 co n 2 and terminate on the 
jp axis ($ * $ * 0) . Figure 20 (a) shows the behavior of a second- 
(order loop (b * 0) with these initial conditions. The terminal or 

jsteady-state phase is sin“^ (0/to n 2 ) as was pointed out in Sec. V. 

! . _ 

Addition of the second integrator causes the steady-state phase to 
become zero with increasing rapidity as b/co n 2 increases. The peak 
phase error also diminishes for increasing values of b/co n 2 at the 
cost, however, of decreased system damping which results in increased 
rms phase error (Figs. 20 (b) - 20 (f)). Finally, as b/co n 2 approaches 

unity the loop becomes unstable. 

The advantages secured by increasing the order of the loop are 

r\ 

more evident in Fig. 21. Here C = -707 again but D/co n =1. It 
was shown in Sec. V that for this and higher values of linear 
frequency shift the loop was incapable of tracking. Figure 21 (a) 
(b/co n 2 * 0) substantiates this fact. On the other hand, even the 
smallest magnitude of th§ second integrator produces zero steady-state 
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phase error. The peak-*phase error is decreased by increasing b/oo n 2 
but as b/co n 2 approaches unity the rms error increases until near 
fc. co n 2 = 1 the system is again unstable. 

✓ 

Similar behavior is apparent in Figs. 22, 23, and 24 except that 
as D/a> n 2 increases an increasingly greater value of b/co n 2 is required 
to keep the loop in lock. Finally, when D/co n 2 is near 2 a value of 
b/<D n 2 near unity is required. However, it is evident from Figs. 20(g) 
to 24 (f) that for this value of b/co n 2 the system is unstable. The 
range of values that b/a> n 2 may take on in order for the loop to remain 
in lock is plotted as a function of D/co n 2 in Fig. 25 - 27 for £ * 0.707, 
0.5, and 1 respectively. The cross-hatched area represents permissible 
values of b/co n 2 . Clearly, for linear frequency shifts the third-order 
loop has extended the tracking range to approximately twice that of 
the second-order loop for a damping factor, C» of the latter equal to 
6.707 and even further for higher values of 

A theoretical explanation of the oscillatory behavior for 
b/»n near or above unity is possible. The characteristic polynomial 
of the loop (Eq. 52) can be rewritten as 


( 2£ co n p + 






n 


T 


<p = o 


(55) 


The Routh-Hurwitz test on the polynomial yields the fact that for 


b sin 

2 ^ 

ton q> 


(56) 


the loop is unstable. This indicates that if in a particular tra- 
jectory the phase error reaches such a value that condition (56) is 
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satisfied the trajectory will become unstable* Clearly, if b/co ft 2 * 1 
-.he condition is satisfied for all 9 and the system is always unstable. 
This fact is also borne out by the linear model which follows from 
the approximation sin <p ca 9. Application of the Routh-Hurwitz test 
to the linear equation sets the instability criterion as b/co n 2 > 1. 

The tracking behavior of the third-order loop does not seem to 
be much different than that of the second-order loop for D/a> n 2 < 1. 

For D/a> n > l, the region for which the second-order loop can not 
lock in, there is a limited range of D/co n 2 for which the proper 
choice of b/co n 2 will secure lock. The pull-in behavior when the VCO 
is not initially tracking the reference signal seems to be less stable 
than for a second-order loop. However, generalizations can not be 
made because of the complexity of the phase-space trajectories. 


VIII. NOISE AND LOOP NOISE ^BANDWIDTH 

Some comments are in order on the effect of noise on the pull-in 
and tracking behavior. A small amount of phase and frequency jitter 
due to noise will not affect the phase-plane trajectories significantly; 
it will result only in hash being superimposed on the trajectories. 
However, the bias due to the signal frequency error is still the 
dominating factor which causes the loop to pull in and a small jitter 
on the trajectories will not affect the number of cycles required to 
achieve lock. On the other hand, if the input noise is significant, 
particularly when the initial frequency error is large, the disturbance 
on the trajectory may be such that a different number of cycles will be 
required for the system to achieve lock. 
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When the loop is in lock and tracking a noisy signal it is 
convenient to consider the linear model which results from the 
approximation sin <p a ?. This enables one to analyze the propagation 
of noise through the system and even to derive the form of the loop 
filter which will result in a Wiener optimum linear system (see Ref. l) 
A key parameter which results from the consideration of noise in the 
system is the loop noise -bandwidth B L defined as 




(57) 


where H(p) is the transfer function of the linearized loop from input 
phase to the output phase of the VCO. The basis of this definition 
lies in the fact that a flat noise spectral density of magnitude <t> 
at the input of the linear system will produce at the output a noise 
power 2«>B l . Thus B L is the bandwidth of the ideal square-cutoff 
low-pass filter which produces the same amount of noise power at its 

output as does the linear system whose realizable transfer function 
is H(p). 

With this definition, the loop noise-bandwidths of the various 
filters nay be considered. For the first-order loop the linearized 
transfer function follows from Eqs. (3) and (4). 


H(p) 


. Tq(p) _ K 

? s (p) P + K 


(58) 


where <p 0 “ ?s ‘ ¥ is the output phase of the VCO and F(p) * l. 
Applying the definition (57) results in a loop noise -bandwidth for 
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the first-order loop of 


Bi 


K 

4 


(59) 


Then the transfer function (58) can be expressed in terns of B^: 


H(p) 


4 B l 
P + 4 B L 


(60) 


Rechtin has shown (Ref. 1) that to minimize the output noise 
power and simultaneously to maintain the rms phase error due to a 
step change in frequency within a given bound , the optimum linear 
filter will be of second order and of the form 



(61) 


From the denominator which is the characteristic polynomial it is 
clear that the damping factor of this optimum loop is C 3 l/*/2* = .707. 
Application of Eq. (57) yields a 


b L 


3o>, 


n 


47 * 


(62) 


Then 


1 + 


H(p) 


4 Bi 


1 + 


4 B L 


32 Bt 2 P 


(63) 
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Also it is shown that to minimize the output noise power and 
simultaneously to maintain the rms phase error due to a linearly 
varying frequency input within a given bound, the optimum filter will 
be of third order and its transfer function in terms of the loop noise 
bandwidth is 


H(p) 


5 25 

1 + 7 — p + ; 

3 Bl 18 B^ 


25 


1 + — - p + 

3 B l 18 B l 


0 2 + 
2 P + 


125 


216 B L ' 


(64) 


The characteristic polynomial of the system is then 


p3 + J B L p2 + S B L 2 P + X25 ®L 3 


Comparison of this with that of Eq. (55) for which the linearizing 
assumption is made results in the following values for the C* co n , an< * 
b parameters in terms of B L . 


co 


n 


6 72 


B 


L 


C 55 * .707 

b = £ b l 2 

25 L 
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(65) 


IX. SUMMARY OF RESULTS 


The various conclusions relative to the filters and signals 
considered are summarized in terms of both the notations which have 
been introduced 

1. First-Order Loop (Sec. III.) 
a. Pull-in range 


Q < K ^ or Q < 
sec 


4 B l 


rad 

sec 


b. Pull-in time for initial frequency error, Q 


,?final 


^ initial 


Q - K sin <j> 


^ final 


^initial 


Q - 4 B'l sin 9 


The integral is evaluated in the standard tables. 

2. Second-Order Loop (Perfect Integrator); Constant Reference 
Frequency (Sec. IV.) 

a. Pull-in range is infinite 

b. Pull-in time for large initial frequency error, Q 
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t etf 


_s!_ 

2 C °>„ 3 


se**- 


A more accurate expression for smaller Q is given 
by Eq. (28) 

For optimum loop (using Eq. 62) 


t P 


27 Q 2 

s sec . 

256 B l 3 


3 . Second-Order Loop (Perfect Integrator); Linearly Varying 

Reference Frequency or Constant VCO Drift at Rate D rad/sec 2 
(Sec. V) 

a. Approximate pull-in range for Q/co n 2 < 1/2 

q < rad 

D set 


for optimum loop when D/B L 


2 


< 


16/9 



Q < 


128 B^ 

27' D 


rad 

sec 


b. Maximum rate which can be tracked 


D 



rad 

2 

sec^ 


for optimum loop 


D 


32 n 2 rad 
-9 B L 2 

sec 
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4. Secnnd-Ordsr Leap ( Imperfect Integrator) Constant Reference 
Frequency (Sec. VI ) 

Pull-in range 

Q < 2 co n Vl^T7rad 
n o sec 


where a is inverse time constant of filter. 
For optimum loop 



5. Third-Order Loop (Sec. VII) 

0 

Tracking range can be extended beyond D/co n 2 * 1 by means of a 
second integrator in the loop filter. Figure 25 indicates that if in 
a second order loop, whose damping factor £ =0.707, a second integrator 
is introduced, the tracking range can be extended to D/o& n 2 = 1.83 by 
making the normalized value of the second integrator constant 
b/con^ = 0.65. 

The optimum third order loop in terms of the other set of 
paraaeeters has C “0.707, © n = 6 -72/5 B L and b/co n 2 * 1/4. Its 
tracking range extends to 



1.5 



at 


75 

144 
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The reason that b/te^ i s chosen in this case to yield a non- 
maximum tracking range is that the optiiauiB loop must strike a 
compromise between tracking accuracy and output noise power, so that 
optimizing 8 composite of the two is preferrable to optimizing either 
one. 
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A - A * saddle point asymptotes 

a = first-order integrator gain 

i 

B - B = periodic limit of lock-in 
B l = loop noise-bandwidth 
j b * second-order integrator gain 

« reference signal amplitude 
Cg = VCO output amplitude 
‘ D = linear rate of frequency variation 

J 

i E ss » steady-state phase error 
e^ * multiplier output 
F(p) = filter transfer function 

1 

H( p ) * loop transfer function 

i ^0 = ^ 1^2 

K, = filter DC gain 

i A 

j K 2 * VOO gain 

K = loop gain 
p * Laplace operator 
t = time variable 
| X = 9 phase error 

y * 9 frequency error 

1 

! y A 31 value of y along A - A 
I yg =5 initial y 

y av = average value of periodic trajectory 

a - inverse time constant of filter with imperfect integrator 

P * linear decrease in y per cycle of x - 

j. leat line of .text- or footnote 

! 

;01ASai?J CATION (j.f &*y) 





C “ loop damping factor 
x * normalized time variable 
9 * phase error 

= reference signal phase 
* VOO output phase 
oo * ^ frequency error 
®S * $s reference signal frequency 
® c * VOO center frequency 
oo^ = loop natural frequency 
= initial VOO frequency 
Q ^ cog - coc 

Qq * initial Q 


I 
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Fig. 4. C * 1.414 (ovardawped) 
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Fig. 12. „ C “ 0.707, D/co^ = 0.95, 

* 72%-* 2/5 *?'rad 
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Fig. 14, C * 0.707, <x/2Cco n * 0.1, 

aQ/co^ = 0.4 
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Fig. 15. £ =* 0.707, a/2Cco n = 0.1, 

aQ/co^ = 0.6, Q/co n = 3-7^2 
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Fig. 16. C * 0.707 , a/2Co n * 0.1, 
aQ/co^ = 0.7, Q/co n = 7/2 J2 
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Fig. 17. C * 0.707, a/2Cco^ * 0.1, 
aQ/co^ = 0.8, Q/co n = 4-/2~ 
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Fig. 18. C * 0.707, o/^a^ « 0.1, 
aQ/co* = 0.9, Q/co n = 9/272 
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